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A condensation transition was predicted for growing technological networks evolving by preferen- 
tial attachment and competing quality of their nodes, as described by the fitness model. When this 
condensation occurs a node acquires a finite fraction of all the links of the network. Earlier studies 
based on steady state degree distribution and on the mapping to Bose-Einstein condensation, were 
able to identify the critical point. Here we characterize the dynamics of condensation and we present 
evidence that below the condensation temperature there is a slow down of the dynamics and that a 
single node {not necessarily the best node in the network) emerges as the winner for very long times. 
The characteristic time t* at which this phenomenon occurs diverges both at the critical point and 
at T — > when new links are attached deterministically to the highest quality node of the network. 
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Condensation phenomena [l|, l^i S Ij, la, IS, |7|] in complex 
networks [SJ, |9|, [lOJ, [ill are structural phase transitions in 
which a node grabs a finite fraction of all the links of 
the network. This phenomenon is of particular interest 
in the case of technological networks. The maps of the 
Internet at the Autonomous System Level show that the 
fraction of nodes connected to the most connected node 
is increasing in time reaching a share the order of 10% in 
recent maps [12| . Also in the World-Wide Web the share 
of webpages linked to Google webpages is of the order of 
1% a large number if one takes into account the size of 
the World- Wide- Web. Are the Internet and the World- 
Wide- Web close to a condensation transition? Which are 
the dynamical signatures of a condensation? Which are 
the consequences of a condensation of technological net- 
works? The problem might have relevant implications 
for the monitoring of technological networks and there 
might be important difference between the statistics of 
lead change below and above [l3| the condensation tran- 
sition. In the following we are studying these problems 
focusing on the dynamics of the fitness model [3, [ij] in 
which nodes acquire links in proportion to their connec- 
tivity and in proportion of their fitness which indicates 
the quality of the node. This model, has been considered 
a good stylized model for the Internet ll| and is a good 
stylized model describing the emergence of high quality 
search engines in the World- Wide- Web and show a con- 
densation phase transition depending on the parameters 
of the model. The study of the dynamical properties of 
this model will be able give some estimates of the char- 
acteristic time scale at which a condensate node might 
emerge in a condensation scenario and will allow us to 
evaluate, once a condensate node is formed, the proba- 
bility that a new node would overcome the condensate at 
later times. 

Questions concerning the dynamics of condensation in 
within the fitness model, might also shed some light on 



the relation between the condensation phenomena occur- 
ring in this model an the other off-equilibrium conden- 
sation phase transitions occurring in complex systems ( 
traffic jam, wealth distribution, urn and network models) 
S \m [13, H E, M, [U, m. Hi . Off-equillbrium phe- 
nomena have attracted large interest in the last ten years. 
Mainly the attention has been addressed to the character- 
ization of the steady state of these out-of equilibrium sys - 
tcms, above the condensation transition [2, y, lZ|, ll6|, ll7| . 
On the cont rary th e dy namics by which the condensate 
emerges IT], llSl, ll9|, \2^ has been studied only in the case 
of urn models and models which can be reduced to the 



Zero Range Process [15|. However the models where the 
condensation occurs as a structural phase transition of a 
network arc not in general reducible to the Zero-Range- 
Proccss. 

In this paper we study the dynamics of condensation 
in the fitness model. We apply the rank statistics for 
the description of the dynamics of condensation and we 
will show that below the condensation phase transition a 
condensate node emerges only after a characteristic time 
t* . This node is then one of the nodes with highest fit- 
ness. The probability that later-comers high fitness nodes 
will overcome the condensate is decaying with time and 
there is a slow down of the dynamics. At the condensa- 
tion transition the characteristic time t* diverges and the 
condensation is only marginal, i.e. we have a sequence of 
highly connected high-fitness nodes each one overcoming 
the other each of then grabbing a very small fraction of 
the total links of the network. 

The fitness model and the condensation phase 
transition- The fitness model |14l | is a growing network 
model. We start from a finite connected network of A^o 
nodes. At each time ti a new node i and m < Nq new 
links are added to the network. To the node i it is as- 
signed a quenched variable, £i ('energy' of the node ) 
drawn from a g{e) distribution. The variable e^ indicates 



the intrinsic quality of the node (lower 'energy' better 
quality). Each of the m new links that are added to the 
network at time <;, connects the new i node to a node j 
chosen with probability 



n,(t) 






(1) 



where we have introduced the fitness r]j associated to 
node j as 



Vj 



-Pej 



(2) 



The parameter /3 = l/T tunes the relevance of the quality 
of the nodes in the choice of the target node. In the limit 
case T — > oo the dynamics of the attachment of new links 
is only dependent on the connectivity of the nodes and we 
recover the Barabasi- Albert model [2J]. For T <C 1 the 
intrinsic quality of the nodes highly affects the dynamic 
of the system. 

For the g{e) distribution such that g{£) ^ as e — > 
a phase transition can occur at a critical temperature 
Tc- Above this critical temperature every node has an 
infinitesimal fraction of all the links, below the critical 
temperature one node grabs a finite fraction of all the 
links. These results can be obtained by consideration 
of the characteristics of the steady state of the model 
above the phase transition. In fact the model can be 
solved in a mean-field approximation by making a self- 
consistent assumption. In particular in 1J| it is assumed 
that Zt = J2j Vjkjit) -^ (Zt) ^ me-f^H + 0{t°). With 
this assumption it can be shown that, at time i, the av- 
erage connectivity of node i arrived in the network at 
time ti grows as a power-law {ki(t)) = m {t/tiY^^^' with 
/(e) = e^^'^^^f^^ and with the constant /i determined by 
the self-consistent equation 



l=I{|3,^i)= deg{e) 



1 



;/3(e-M) _ 1' 



(3) 



In [2] it was shown that the self-consistent equation ([3|) 
is equivalent to the equation for the conservation of the 
total number of links present in the network 



2mt = mt + mtl{j3, /i). 



(4) 



If Eq. ([3]) has a solution then the degree distribution of 
nodes of fitness 77 = e~^'^ is given by Pe{k) ~ fc-i//(-)^i 
and P{k) = J deg{e)Pe{k) is dominated by the term 
Po(e) = k~'' with 7 = e~^'^ + 1 but can have logarithmic 
corrections [IJ]. Nevertheless, if /(/3, 0) < 1 the self- 
consistent equation ([3]) cannot be solved. In this case 
the self-consistent approach fails and we have the con- 
densation phase transition. A necessary condition for 
condensation, ( i.e for /(/3,0) < 1) is that the distribution 
g{e) — » for e ^ 0. For this type of g{e) distributions 
there can be a critical temperature Tc = l//3c such that 



I{/3, 0) < 1 V f3 > (3c- This phase transition is for- 
mally equivalent to the Bose-Einstein condensation for 
non-interacting Bose gases in dimensions d > 2. Using 
the similarity to the Bose-Einstein condensation in Bose 
gases, in Ref. [2[ it was assumed that below the conden- 
sation temperature the equation for conservation of the 
links can be written as 



2mt = mt + mtliP, 0) + nQ{f3)mt 



(5) 



where no{P) indicates the number of links attached to the 
most connected node. Consistently with ([5|), simulation 
results reported in [3| show that indeed the fitness model 
undergoes a condensation phase transition at T = Tc- An 
example of (7(e) distributions for which there is a conden- 
sation is 



g{e) = {0 + l)e' 



(6) 



and £ € [0, 1]. Here and in the following we will always 
consider g{e) distributed according to dH). Assuming ([5]) 
and a distribution of the type ^ the fraction of con- 
densed links, similarly to the corresponding result for the 
Bose gas, will take the form 



noif3) 



e+i 



(7) 



Dynamics of the fitness model- The necessary condition 
for condensation g{e) — > as e — *■ is a clear indica- 
tion that dynamical effects not captured by the station- 
ary state approaches [3, y, 0] will be very important for 
the emergence of a long lasting condensate node in the 
network. Moreover all considerations regarding the dy- 
namical aspects of the condensation phenomena on the 
fitness network cannot take advantage to the similarity 
with the condensation phenomena in the non-interacting 
Bose gas. In fact, in a Bose gas, the equivalent of no(/3) 
is the occupation of the state at energy e = 0. On the 
contrary, in the fitness networks, no{(3) cannot be the 
connectivity of the fittest node (e = 0). Indeed, in the 
fitness model, the fittest node appears in the network 
only in the infinite time limit because the probability 
g{e) of having a node of quality e, goes to zero as e — > 0. 
Therefore the question is: is the condensate node the best 
node in the network? Once the condensate has appeared 
in the network is there a possibility for later-comers high 
fitness nodes to take over the condensate? To narrow 
down this general question we define a specific problem 
we want to address. First of all, let us introduce the no- 
tion of a record. A record en{t) at time t is the node 
(or equivalently the value of its e) with minimal e (max- 
imal fitness) present in the network. In the fitness model 
below the condensation transition, all records are good 
candidates as condensate nodes. Suppose that a record 
e, arrived in the network at time i, is a condensate with 
k links at the time tr where a the r-th record occurs. 




FIG. 1: We plot the connectivities of records as a function 
of time. The time tr indicates the time at which a record 
node with Sr arrive in the network. The effective model only 
focuses on the competition for links of the most connected 
nodes with the node of highest fitness. 



Wc want to know the probability that the r-th record of 
fitness r]r = e~^^^ overcome the condensate node before 
the subsequent record e^+i enters the network (See Fig. 
[T] for a pictorial representation of our problem) . In order 
to answer this question we have to use some results re- 
lated to record statistics. 
Records- The typical value of the record at time t is given 

by 



1 

g{x)dx = -, 



(8) 



indicating that the average number of nodes with e > 
£r{T) in a network of t nodes is less than one. Using ([5]) 
for the distribution ([6| we obtain 



efl(t) ~< "+1. 



(9) 



The statistics of records is a field on its own [25|, here 
we cite only another result which will be used in the fol- 
lowing: given the r-th record at time tr the probability 
to have the following record at time f,._)-i > t in the ap- 
proximation of a continuous time is given by [2a | 



PRiU+l >t,tr) = j 



where t > t,.. 

The effective model- At low temperatures, below the con- 
densation phase transition, the linking probability IIj 
(Eq.([T])) of the fitness model will be relevant only for 
nodes of very high fitness (the series of records) or for the 
condensate node. Therefore we consider the extremely 
simplified effective model of just two competing nodes 
(the condensate node and the record) disregarding all 
the other nodes present in the network. This simplifica- 
tion will provide us a scenario for condensation which we 
will subsequently compare with simulation results on the 
original fitness model. At time tr there is a node (the 
condensate) arrived in the network at time i with fitness 



77 = e~^^ and k links, and a record node with fitness 
r]r = e~^^'' which is just born and has degree 1. At each 
time we distribute one link to one of the two nodes with 
probability given by expression ([T]). 

The probability pt(r) that the node with fitness rjr has 
T links at time t satisfies in the continuous time limit 



dt 



VrJT-l) , , IJrT , , , . 

-— Pt{T-l)-^Pt{T), (11) 



where the competition is included in the formula through 
Zt^T = vik -\- t — tr — t) -\- rjrT. We are interested in 
the dynamics of this effective model until the node with 
fitness rjr takes over, i.e. until 



T <Tc ^ {k + t- tr)/2. 



(12) 



Equation (jlip is non-trivial to study and some approxi- 
mations are necessary. Therefore we assume that Zt^r — 
Zt,T-i, (i-e. f3e <C 1) which, taking into account the scal- 
ing of the records with the time of their appearance, Eq. 
(IHl) reduces to the condition 



i->r-(9+l) 



(13) 



Furthermore we assume that C ~ Zt^r/(k + t — tr) 
can be approximated by a constant for t < tc- In- 
deed for T < Tc we have that C varies in a narrow 
interval C G {fj, {fj + rjr)/ 2). In this approximation we 
can rewrite equation (jllll using the generating function 
q{z,t) = X]t=iP(''''^)'^^- ^^'^ equation becomes 

CCk + t-tr)^^^VrZ{z-l)^^ (14) 

which, if we impose the initial condition q[z, tr) — z, has 
the solution 



q{z,t) 



( k+t- 






-r,r/C 



+ (1-^) 



k + t~tr 



-Vr/C 



and the expression for the probability pt (r) becomes 



(10) pt{r) 



rZ "7" t ^7" 



flr/C 



k + t-tr 



vjc\ 



The probability that, at some time t, t < tq, i.e. the 
probability that the record has not become the most con- 
nected node of the network can be calculated: 



Pi {t <Tc) ~ 1 - exp 



where in the last step we have assumed that tc — tr ^ 1 . 
Expression P^ for Pi(r < tc) allow us to evaluate the 
characteristic time tc at which the record r becomes the 
most connected node, i.e. 
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iv,-c)/c- 
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\k-\-t-tr, 



tc-tr~¥+^/^^^-^\ 



(15) 
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FIG. 2: Probability Pc that a record overcomes the most 
connected node of the network in the fitness model averaged 
over 10'' runs for a time lapse of tmax = 10^ . 



FIG. 3: Characteristic time t* calculated with the estima- 
tion (|20p (B=10) after which a long lasting condensate node 
emerges for the fitness model for 6 = 1 and a = 1, oo 



Therefore, by making use of (|TO|l we can evaluate the 

re 

probabiUty P™ that the record r has become the most 
connected node of the network before the appearance of 
the record r + 1, 



P 



eff 



PFt{tr+l > tc,tr 



nkC/irir^C) 



16) 



where h — k/tr- Taking into account the interval of 
definition of C, the scaling exponent a = C /{rj^. — C) £ 
[l/[f3{e - er)],2/[P{e - Er)]]- Therefore if e > e^ then 



1 



ITi^+T. 



(17) 



In other words as long as Eq. (|T3|) is satisfied, the expo- 
nent a is greater than one, (i.e. a > 1) and it increases 
with larger t times and for higher temperatures T. 

If we fix the condition P^ < 1/B we find a character- 
istic time for the appearance of a long lasting condensate 
node in the effective model is given by 



t* - max T-C'+i) 



B\^/" l' 



(18) 



where n ~ k/t^. 

Dynamics of condensation- We expect that the effective 
model, describing the competition of only two nodes, pro- 
vides an upper bound for the probability Pc{k) that 
in the fitness model a new record overcomes the most 
connected node before the appearance of a subsequent 
record. We define Pc{k) the probability that a record, 
introduced in a network with a condensate of connectiv- 
ity fc, is able to overcome the connectivity of the conden- 
sate node before the subsequent record. To numerically 
evaluate Pcik) we consider all the records that follow 
the emergence of a condensate and we evaluate the ra- 
tio between the number of positive events and the total 
number of records that appear in the network when the 
condensate has k links. The data reported in Figure [2] 



are taken for 6* = 1 at different temperatures T below the 
phase transition. The results shown are statistics taken 
over 10"* networks for a time lapse T~*^^+^^ < t < 10^. 
Figure [2] shows clearly that the probability Pc(k) decays 
as 



Pc{k) ex k 



(19) 



with a > 1 and that the exponent a increases with T 
as predicted by the effective model. Using the effective 
model estimate (|18p for the characteristic time t* for hav- 
ing a long lived winner node which condenses, and assum- 
ing fi ~ no(/3) with ri,o(/3) given by ([7]) we get 



max<^T-(»+i),Bi/" 



-(l + l/o 



(20) 

This characteristic time is depicted in Figure [3] for dif- 
ferent values of a and the results are consistent with a 
divergence of i* at T ^ Tc and at T — > 0. In order to 
evaluate the time life of a condensate for times t > t* we 



can use the scaling of the probability Pc 



'k- 



and 

we can assume that for the condensate node k ~ n^t. 
Moreover if the relevant competition is only between the 
condensate node and the node with highest fitness, the 
probability to have the same condensate after S records 
is given by 



Pf{S) 



r+S 

n 

s— r* + l 



1- 



1 



m{ni)ts 



where ts are the times of record occurrence after the 
condensate has become the most connected node in the 
network (i.e. for s > r*). The average of this quan- 
tity over the distributions of the times of the records 
UsPRits,ts-i)=Usts-i/tl gives 



P#)(5)^cxp 



1 



1 



(1 



1 



ano{notr*)" 



(21) 



Therefore for long times, we expect (P^)(5) -^ const 
as 5* ^ cxD. Consequently we can conclude that for times 
t > t* the condensate node is typically not the best node 
of the network nevertheless it dominates the network for 
very long times. Only for T —> Tc and T —> the time 
needed to havelong lasting condensate diverges and in 
the network there is a fair competition and a succession 
of high fitness nodes on which condensation occurs. 

Conclusions- In conclusion we have studied the dynam- 
ics of condensation in growing network models within the 
fitness model. We show that below the condensed phase 
transition there is a characteristic time t* . For times 
t < t* the network dynamics is dominated by a succes- 
sion of high fitness nodes with a finite fraction of the 
links, above the characteristic time t* a long-lasting win- 
ner takes over acquiring a finite fraction of all the links 
and slowing down the dynamics. The characteristic time 
t* is diverging at the critical point of the condensation 
phase transition and in the limit T — > 0. These results 
may have strong implications for the monitoring of tech- 
nological networks in which we expect a fair competition 
lasting for very long times only if close to the condensa- 
tion transition. 
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